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We investigate how the frequencies of gravity modes depend on the detailed properties of the chemical composition 
gradient that develops near the core of main-sequence stars and, therefore, on the transport processes that are able to 
modify the p profile in the central regions. We show that in main-sequence models, similarly to the case of white dwarfs, 
the periods of high-order gravity modes are accurately described by a uniform period spacing superposed to an oscillatory 
component. The periodicity and amplitude of such component are related, respectively, to the location and sharpness of 
the p gradient. We briefly discuss and interpret, by means of this simple approximation, the effect of turbulent mixing 
near the core on the periods of both high-order and low-order g modes, as well as of modes of mixed pressure-gravity 
character. 
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1 Introduction 

In main sequence stars the periods of gravity modes are sen- 
sitive probes of stellar cores and, in particular, of the chem- 
ical composition gradient that develops near the outer edge 
of the convective core. 

High-order gravity modes are observed in two classes 
of mai n-sequence star s: 7 Doradus (M ~ 1.4 — 1.8 M Q , 
see e.g. iHandlerll 19991 for a review) and Slowly Pulsating B 
stars, with masses from about 3 to 8 Mq and spectral type 



B3-B8 (rWaelkenslll991l) . The seismic modelling of these 
classes of pulsators is a formidable task to undertake. The 
frequencies of high-order g modes are in fact closely spaced 
and can be severely perturbed by the effects of ro tation (see 
e.g. lDintrans & Rieutordl200filSuarez et al.l2005l) . Nonethe-| 
less, the high scientific interest of these classes of pulsators 
has driven efforts in both the observational and theoreti- 



interaction between p and g modes. In Sec.[3]we give an ex- 
ample on how different mixing processes acting near the 
core affect differently the oscillation spectra observed in 
these stars. 



2 High-order g modes 

As it is well known, the period spectrum of gravity modes is 
determined by the spatial distribution of the Brunt-Vaisala 
frequency (N) which is defined as: 
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cal domains. As suggested by ISuarez et al.l (I2005h in the 
case of 7 Dor stars, a seismic analysis becomes feasible 
for slowly rotating targets. In the se favorable c ases the first- 
order asymptotic approximation ( Tassoulfl980l) can be used -/V 2 — 
as a to ol to derive the buoyancy radius of the star (see Moya et aH 
2005) from the observed frequencies. The g-mode spectra of where 



where p is density, p is pressure, g is the gravitational accel- 
eration and Ti = (dlnp/dlnp) d . 

N can be approximated, assuming the ideal gas law for 
a fully-ionized gas, as: 
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these stars contain however much more information on the 
internal structure of the star: this is illustrated in Sec.|2]by 
means of a simple refi nement of the asymptotic approxima- 
tion of lTassoull(ll980l) . 

In Sec. [3] we briefly discuss the case of low-order g 
modes, and modes of mixed pressure and gravity character. 
In main-sequence stars these modes are dete cted in S Scuti 
(M 2 ± 1.5—2.5 M Q , see e.g. the overview by Rodriguez & Brei 
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) and Cephei (M ~ 7 - 20 M , Stankov & Handler 



2005) pulsators. Moreover solar-like oscillations in subgiant 



stars can also present a p-mode spectrum perturbed by the 



T is temperature and p is the mean molecular weight. 

The term V M gives the explicit contribution of a change 
jrf^hemical composition to N. The first or der asymptotic 
approximation developed by Tassoul (1980) shows that, in 
the case of a model that consists of an inner convective 
core and an outer radiative envelope (we refer to the work 
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by iTassoul 1980, for a complete analysis of other possi- 
ble cases), the periods of low-degree, high-order g modes 
(where iV 2 3> u 2 ) are given by: 



Pk = 



Lf 1 

Jx 



N\ 



dx 



(2k + n e ) 



(3) 



where L = [£(£ + l)] 1 / 2 (with i the mode degree), n e the 
effective polytropic index of the surface layer, x the nor- 
malized radius and xq corresponds to the boundary of the 
convective core. In order to avoid confusion with n e , the 
radial order of g modes is represented by k. 

Following Eq. [3] the periods are asymptotically equally 
spaced in k and the spacing decreases with increasing L. It 
is therefore natural to introduce, in analogy to the large fre- 
quency separation of p modes, the period spacing of gravity 
modes, defined as: 




0.2 
m/M 



AP = P k+1 - P k 



(4) 



In this work we show that deviations from a constant AP 
contain information on the chemical composition gradient 
left by a convective core evolving on the main sequence. 

The signatures of chemical stratifications in high-order 
gravity modes have been extensively investigated theoret- 
ical ly, and then ob served, in pulsating white dwarfs (see 
e.g. Kawalerfl995 , for a review). The influence of chemical 
composition gradients on g modes in main sequence stars 
has been partly addressed and suggest ed i n two farsighted 
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Fig. 1 Behaviour of the Brunt- Vaisala frequency (upper left 
panel), of the hydrogen abundance profile (upper right panel) and 



works bv lBerthomieu & Provostld 19881) and lDziembowski et alftf the 1 = 1 g-mode period spacing in models of 6 M . We 

consider several models along the main sequence with decreas- 
ing central hydrogen abundance (X c 0.7, 0.5, 0.3 and 0.1). In the 
lower panel the uppermost model has X c = 0.7 and the lowermost 



( 1993). Inspire d by these works and following the approach 
of iBerfhomieu & ProvosJ dl988l) and lBrassard et al.1 d 19921) 
we investigate the properties of high-order, low-degree grav-| 
ity modes in main sequence stellar models. 

Whil e a detailed descripti on of our investigation is pre- 
sented in lMiglio et al.1 d2007bl) . we here recall the results of 
a simple approximation of g-mode periods that allows to 
relate the departure from the c onstant period spacing pre- 
dicted by the approximation by lTassoull dl980l) . to the prop- 
erties of the variation of the Brunt-Vaisala frequency lo- 
cated in the /i-gradient region near the core. 

This simple approach makes use of the varia tional prin- 



X c = 0.1. 



imation, that reduces the differential equations of stellar adi- 
abatic oscillations to a system of the second order. Further- 
more, since we deal with high-order gravity modes, the eigen-| 
function s are well des cribed by their JWKB approximation 
(see e.g. lGoughlll993l) . 

We then model the sharp feature in ^ located at x = 
X,, as: 



ciple for adiabatic stellar oscillations (see e.g. lUnno et al 
Il989l) . The effect of a sharp feature in the model (a chemical 
composition gradient, for instance) can be estimated from 
the periodic signature in SP, defined as the difference be- 
tween the periods of the star showing such a sharp variation 
and the periods of an otherwise fictitious smooth model. 

We consider a model with a radiative envelope and a 
convective core whose boundary is located at a normalized 
radius xq. N- and N + are the values of the Brunt-Vaisala 
frequency at the outer and inner border of the /x-gradient 
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N a 2 
where H(x) is the step function. 

Retaining only periodic terms in SP and int egrating by 



parts w e obtain, after simple calculations (see Migli o et al. 
2007bh : 
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with N+ < 



N-. Then 
0, a sharp 



region. ^....^ « — , N 
a = 1 describes the smooth model and a 
discontinuity in 7Y. 

To obtai n an estimate of SP, we adopt (following the 
approach by Montgomery et al . 2003) the Cowling approx- 



where is a phase constant; the buoyancy radius of the star 
is defined as: 
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and the buoyancy radius of the discontinuity is: 
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Fig. 2 Behaviour of the Brunt-Vaisala frequency (upper left 
panel), of the hydrogen abundance profile (upper right panel) 
and of the 1 = 1 g-mode period spacing in models of 6.0 
Mq computed with a turbulent diffusion coefficient Dt = 5000 



For small 8P we can substitute in the expression above 
th e asymptotic ap proximation for g-mode periods derived 
bv .Tassoull(ll9801) : 



P fc =7T 2 ^(2fc+0') , 

where </>' is a phase constant, and find: 

+ . (9) 

From this simple approach we derive that the signature of a 
sharp feature in the Brunt-Vaisala frequency is a sinusoidal 
component in the periods of oscillations, and thus in the pe- 
riod spacing, with a periodicity in terms of the radial order 
k given by 

Afc ~ 5*i • (10) 
Ho 

The amplitude of this sinusoidal component is proportional 
to the sharpness of the variation in N and does not depend 
on the order of the mode k. 

Such a simple approach allows us to easily test the ef- 
fect of having a less sharp "glitch" in the Brunt-Vaisala fre- 
quency. We model SN as a ramp function instead of a step 
function: 



SN 
~W 
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H(x ll - x) . 



(11) 



In this case integration by parts leads to a sinusoidal 
component in SPk whose amplitude is modulated by a fac- 
tor 1/Pk and therefore decreases with increasing k, i.e. 
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The information contained in the amplitude of the sinusoidal^ 
component is potentially very interesting. It reflects the dif- 
ferent characteristics of the chemical composition gradient 
resulting, for example, from a different treatment of the mix-| 
ing process in convective cores and from considering extra- 
mixing processes acting near the core. 

With this simple approximation we can explain the be- 
haviour of the period spacing in main-sequence stars as a 
superposition of the constant term derived from Eq. [3] and 
periodic components related to V M . 

As an example we consider in FigureQ]the high-order g- 
mode spectrum of 6 M Q models along the main sequence. 
The amplitude and periodicity of the components in AP 
can be easily related to the profile of the Brunt-Vaisala fre- 
quency by means of expression ( fTOb . For instance, the AA: = 
7 in Fig. Q] (lower panel, X c = 0.50), corresponds to the 
sharp signal in N at IIo/II r ~ 0.14. As the star evolves, 
the sharp feature in N is shifted at higher Ho /IT , and when 
Iln/II^ ~ 0.5 a sort of beating appears in the period spac- 
ing due to the fact that the sampling frequency is about half 
the frequency of the periodic component. 

The behaviour of the period spacing is therefore mainly 
determined by the location and sharpness of the glitch in 
in the /i-gradient region. In m a in-sequ ence stars the lat- 
ter, as described in Migl io et al. d2007bh . depends on the 
size of the convective core and on extra-mixing processes 
(e.g. overshooting, diffusion, turbulent mixing) that may al- 
ter the chemical composition profile in the central region of 
the star. 

As an example we consider models of a 6 M star com- 
puted including as extra mixing turbulent diffusion near the 
core (described by a turbulent-diffusion coefficient Dj- = 
5 x 10 3 cm 2 s _1 , see Montalban et al., these proceedings 
for a description of this simple parameterization). In Fig. [2] 
we notice that such a mixing has a substantial effect on the 
period spacing: the amplitude of the periodic components 
in AP becomes a decreasing function of the radial order 
k. This behaviour can be easily explained by the analytical 
approximation presented (see Eq. [T2l . provided the sharp 
feature in N is modelled not as a step function but, for in- 
stance, as a ramp. 

3 Low-order g modes and mixed modes 

Among the different classes of main-sequence pulsating stars,| 
low order gravity modes and mixed-modes are found to be 
excited in j3 Cephei and 8 Scuti stars. Moreover, in the case 
of stars evolving as sub-giants, mixed modes are also ex- 
pected in the frequency domain typical of solar-like oscilla- 
tions. 
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1.5x10" 



Fig. 3 Behaviour of the hydrogen abundance profile (upper 
panel) and of the Brunt- Vaisala frequency (lower panel) in 10 
Mq models with X c ~ 0.3. The different lines correspond to 
models calculated with overshooting (qov = 0.2) (dashed) and 
turbulent mixing with Dt = 2.8 x 10 5 cm 2 s _1 (continuous). 



In the previous section we recalled that the properties of 
high-order g-mode spectra can be easily related by means 
of an analytical approximation to the detailed characteris- 
tics of the "glitch" of TV located in the /i gradient region. 
This approximation, however, becomes invalid when con- 
sidering modes of low-order, and cannot be used to give an 
accurate description of the properties of such modes. A fur- 
ther complication to an analytical description of low-order 
non-radial modes is given by the interaction between pres- 
sure and gravity modes that gives rise t o modes of mixe d 
pressure and gravity character (see e.g. Shibahashi 19791) . 
Nonetheless, as we describe in the following, the analytical 
description of gravity modes presented in Sec.|2]is still able 
to qualitatively relate some properties of g-mode spectra for 
modes of low-order, to the characteristics of the /i-gradient 
region near the core. 



3.1 /3 Cephei and S Scuti stars 

In main-sequence stars, the frequency of a mixed mode is 
known to provide information about the extension of the 
central mixed region. It has been shown, for instance, that 
the detection of the frequency of a mixed mode would al- 
low the determinati on of the overshooting parameter in 5 
Scuti stars (see e.g IPamvatnvkhll 19991: iMichel et all 1 1993 




Audard et al.ll 19951) - A similar inference is also possible in 
j3 Cep stars, where the observed oscillation frequencies are 
now used to place constraints on extra-mixing at the border 

of the convective core (see e.g . Aerts et al.l2003l:|Pamvatnvkh et^ af.jf 1 ^ var i at i° n m N 
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Fig. 4 Period spacing for t = 1, 2 and 3 low order modes (up- 
per, central and lower panel) as a function of the period for 10 
Mq models with X c ~ 0.3 computed with overshooting qov = 
0.2 (triangles connected with solid lines) and with turbulent mix- 
ing Dt = 2.8 x 10 5 cm 2 s -1 (squares connected with dashed- 
dotted lines). Horizontal short-dashed lines represent the asymp- 
totic value of the period spacing for gravity modes ( AP oc const), 
whereas the long dashed line shows the asymptotic value of AP 
valid for high order pressure modes (AP oc P 2 ). The frequency 
range typical of a Cep star is delimited by vertical dotted lines. 



Sec. |2] we compute models including turbulent mixing near 
the core, which leads to a smoother /i gradient (and thus 
a smoother N, see Fig. [3]) near the boundary of the con- 
vective core. We compare the effects of such a mixing with 
those of overshooting. The models used in the comparison 
are computed with an overshooting parameter chosen such 
as to lead to an evolutionary track on the HR diagram close 
to that obtained with turbulent mixing. 

We here compare 10 M Q models with X c ~ 0.3 com- 
puted with overshooting (aov = 0.2) and with turbulent 
mixing (Dt = 2.8 x 10 5 cm 2 s _1 ). In Fig.|4]we show that 
in the typical j3 Cep frequency domain (vertical lines) the 
difference between the frequencies of these 10 M models 
can still be qualitatively related to the different behaviour of 



2004; MazumdaretaL||2006; BriquetetalJ|2007). 



Here we investigate, as in the case of high-order g modes,| 
whether the oscillation frequencies typical of f3 Cephei and 
5 Scuti stars are sensitive to the detailed shape of the /i gra- 
dient at the edge of the core (and thus to the mixing process 
that generates such a chemical composition profile). As in 



As shown in the previous section, turbulent mixing leads 
to a reduction in the amplitude of the periodic components 
in the period spacing and such a reduction increases with 
the order of the mode (see Fig. [2j. 

In Fig. [4] we compare the period spacing for ^ = 1,2 and 
3 modes of a 10 M Q model computed with turbulent mixing 



(2006) 
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Fig. 5 HR diagram showing evolutionary tracks of 2 
Mq models calculated without extra-mixing processes (dotted 
line), with overshooting (qqv = 0.2, dashed line) and with tur- 



bulent mixing (Dt = 800 cm 



, full line) 



or with overshooting. The decrease of the amplitude of the 
periodic components in AP, due to the effect of turbulent 
mixing, significantly affects the periods of low-order gravity 
modes and of mixed modes as well. We also compare in 
Fig. H]the period spacing of low-order modes and of mixed 
modes to the expression of AP valid in the two asymptotic 
regimes: for high-order gravity modes (AP oc const) and 
for high-order pressure modes (AP oc P 2 ). 

We refer to Montalban et al., these proceedings, for a 
more detailed description of this effect and for a quantitative 
estimation of the the turbulent mixing induced by rotation in 
models of typical /3 Cep stars. 

Though we concentrated our investigations primarily on 
the effects of turbulent mixing on models of (3 Cephei stars, 
we here recall that a similar effect is also expected in the 
case o f pulsation modes of 5 Scuti stars. In Goupil & Talon 
(12002b it was in fact already suggested that "[. . . ] extra mix- 
ing due to rotation modifies the shape of the inner maxi- 
mum of N. Signatures of these fine structure in frequencies 
of mixed and g modes should be detectable." . As an exam- 
ple of this effect we compare models of 2 M Q stars com- 
puted with turbulent mixing (D T = 800 cm 2 s^ 1 ) and with 
overshooting («ov = 0.2). The corresponding evolutionary 
tracks in the HR diagram are shown in Fig. [5] whereas the 
behaviour of the oscillation frequencies as a function of T e g 
is presented in Fig. [6] Like in the case of (3 Cep models we 
notice that turbulent mixing, despite having an effect on the 
evolutionary tracks similar to that of overshooting, affects 
differently the properties of pulsation modes. For a quantita- 
tive investigation of this effect, however, comparisons with 
models where turbulent mixing results, for instance, from a 
consistent treatment of rotation are still needed. 
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Fig. 6 Frequencies of I — 2 pulsation modes as a function 
of log T e g for main-sequence models of a 2 Mq star. Gray dots 
represent the frequencies of models computed with overshooting 
(aov = 0.2), whereas black circles are the frequencies of models 
computed with Dt = 800 cm 2 s" 1 . Excited modes are repre- 
sented by larger dots (models with overshooting) and thicker cir- 
cles (models with turbulent mixing). 



3.2 Mixed modes in solar-like oscillations 



As discussed e.g. in Di Mauro et al. (2004) and Miglio et al. 



(2007a), the detection of mixed modes in the spectrum of 
stars showing solar-like oscillations would allow to deter- 
mine the evolutionary state of the star and, therefore, to dis- 
criminate among theoretical models computed with differ- 
ent values of the overshooting parameter. In fact, in stars 
that are still in the core-hydrogen burning phase, the fre- 
quencies of gravity modes are not high enough to enter the 
domain of solar-like oscillations, whereas this is no longer 
the case in more evolved models, where mixed modes can 
severely perturb the solar-like oscillation spectrum. 

As in the case of more massive stars presented above, 
we expect that the frequency of a mixed mode depends on 
the detailed properties of the chemical composition gradient 
near the core. 

We here consider, as an example, two models of 1.4 
M subgiant stars with the same surface properties (L,T e s) 
but that are computed including different extra mixing pro- 
cesses in the core: overshooting (aov=0-15) and turbulent 
mixing (Dt = 40 cm 2 s _1 ). Even if the models have al- 
ready left the main sequence, we see in Fig|7]that their chem-| 
ical composition profile (upper panel) still bears the signa- 
ture of a different mixing process in the core, which also 
affects the Brunt-Vaisala frequency (lower panel). In anal- 
ogy to Fig.HJ we show in Fig. [8] that the differences in the 
frequencies of gravity modes, due to a different behaviour of 
N, can also affect the spectrum in the domain of solar-like 
oscillations (vertical dotted lines). A more detailed discus- 
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Fig. 7 Hydrogen profile (upper panel) and Brunt-Vaisala fre- 
quency (lower panel) in models computed with overshooting (con- 
tinuous line) and with turbulent mixing (long-dashed line) 
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Fig. 8 Period spacing for I = 1 modes as a function of the 
period for a 1 .4 M© model computed with overshooting («ov = 
0.15, solid lines) and with turbulent mixing (Dt = 40 cm 2 s -1 , 
dashed-dotted lines). The long-dashed line represents the asymp- 
totic value of AP valid for high order pressure modes. The fre- 
quency range of solar-like oscillations is delimited by vertical dot- 
ted lines. 



sion on this effect and, in particular, on its magnitude and 
detectabil ity for the specific case of the binary 12 Bootis is 
given in lMiglio etail J2007ah . 



4 Conclusions 

We have shown that the periods of gravity modes in main- 
sequence stars can be related, by means of analytical ap- 
proximations, to the detailed characteristics of the /i-gradienfl 
region that develops near the energy generating core and 
thus to the mixing processes that affect the behaviour of /i 
in the central regions. 



Further investigations are however needed to assess un- 
der which observational conditions such information can be 
recovered from the oscillation frequencies, given realistic 
observational errors and, in the case of SPB and 7 Dor pul- 
sators, given the severe influence of rotation on the spectrum 
of g modes. 
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